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Abstract 

A new quantum mechanical notion — Conditional Density Matrix — proposed by 
the authors ||5|, is discussed and is applied to describe some physical processes. This 
notion is a natural generalization of von Neumann density matrix for such processes 
as divisions of quantum systems into subsystems and reunifications of subsystems into 
new joint systems. Conditional Density Matrix assigns a quantum state to a subsystem 
of a composite system under condition that another part of the composite system is in 
some pure state. 

1 Introduction 

A problem of a correct quantum mechanical description of divisions of quantum systems into 
subsystems and reunifications of subsystems into new joint systems attracts a great interest 
due to the present development of quantum communication. 

Although the theory of such processes finds room in the general scheme of quantum 
mechanics proposed by von Neumann in 1927 even now they are often described in a 
fictitious manner. For example, the authors of classical photon teleportation experiment p[ 
write 

The entangled state contains no information on the individual particles; it only indicates 
that two particles will be in the opposite states. The important property of an entangled pair 
that as soon as a measurement on one particles projects it, say, onto 

I <->■> the state of the other one is determined to be \ |>, and vice versa. How could a 
measurement on one of the particles instantaneously influence the state of the other particle, 
which can be arbitrary far away? Einstein, among many other distinguished physicists , could 
simply not accept this "spooky action at a distance". But this property of entangled states 
has been demonstrated by numerous experiments. 



2 The General Scheme of Quantum Mechanics 

It was W.Heisenberg who in 1925 formulated a kinematic postulate of quantum mechanics 
[||. He proposed that there exists a connection between matrices and physical variables: 

variable T matrix {F)mn- 

In the modern language the kinematic postulate looks like: 

Each dynamical variable of a system S corresponds to a linear operator F in Hilbert 
space 7i 

dynamical variable T <^==^ linear operator F. 

The dynamics is given by the famous Heisenberg's equations formulated in terms of 
commutators. 

To compare predictions of the theory with experimental data it was necessary to under- 
stand how one can determine the values of dynamical variables in the given state. W.Heisenberg 
gave a partial answer to this problem: 

If matrix that corresponds to the dynamical variable is diagonal, then its diagonal elements 
define possible values for the dynamical variable, i.e. its spectrum. 

{F)mn = fmSmn {fm} is SpCCtrum T. 

The general solution of the problem was given by von Neumann in 1927. He proposed 
the following procedure for calculation of average values of physical variables: 

<jr> = Tr{Fp). 

Here operator p satisfies three conditions: 

1) P+ = P, 
2) Trp = 1, 

3) yip en < ip\pip > > 0. 

By the formula for average values von Neumann found out the correspondence between 
linear operators p and states of quantum systems: 

state of a system p <^==^ linear operator p. 

In this way, the formula for average values becomes quantum mechanical definition of 
the notion " a state of a system" . The operator p is called Density Matrix. 
From the relation 

(< jr = Tr{F+p) 

one can conclude that Hermitian-conjugate operators correspond to complex-conjugate vari- 
ables and Hermitian operators correspond to real variables. 

F ^ J^* ^ F+, 
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jr = JF* <^ F = F+. 

The real variables are called observables. 

From the properties of density matrix and the definition of positively definite operators: 

F+^F, y^/jen <i}\F'il)> >0, 

it follows that the average value of nonnegative variable is nonnegative. Moreover, the 
average value of nonnegative variable is equal to zero if and only if this variable equals zero. 
Now it is easy to give the following definition: 

variable T has a definite value in the state p if and only if its dispersion in the state p is 
equal to zero. 

In accordance to general definition of the dispersion of an arbitrary variable 

V{A) = <A^> - {<A>f, 
the expression for dispersion of a quantum variable T in the state p has the form: 

where Q is an operator: 

Q = F-<T>E. 

If !F is observable then is a positive definite variable. It follows that the dispersion of T 
is nonnegative. And all this makes clear the above-given definition. 

Since density matrix is a positive definite operator and its trace equals 1, we see that its 
spectrum is pure discrete and it can be written in the form 

P = ^PnPn, 

n 

where P„ is a complete set of self-conjugate projective operators: 

Pfi Pm Pm-Pn ^mnPmy ^ ^ Pn P'- 

n 

Numbers {pn} satisfy the condition 

Pn =Pn, < Pn, J^Pn P^Pn = 1- 

n 

It follows that p acts according to the formula 

= J2Pn E ^na{(l>na\^)- 

The vectors (/)na form an orthonormal basis in the space 7i. Sets A„ = {1, A;„} are defined 
by degeneration multiplicities kn of eigenvalues Pn- 

Now the dispersion of the observable T in the state p is given by the equation 



^p(-^) = ^Pn ll'50na| 
n aGAn 
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All terms in this sum are nonnegative. Hence, if the dispersion is equal to zero, then 

if Pn 7^ 0, then Q(t)na = 0. 

Using the definition of the operator Q, we obtain 

if Pny^O, then F(f)na ^ <f>na{F). 

In other words, if an observable T has a definite value in the given state p, then this 
value is equal to one of the eigenvalues of the operator F . 
In this case we have 

pF(t>noc = (t>naPn{^) , 
Fp4>na = Ka{^)Pn, 

that proves the commutativity of operators F and p. 

It is well known, that if A and B are commutative self-conjugate operators, then there ex- 
ists self-conjugate operator T with non-degenerate spectrum such that A and B are functions 
of f: 

na 

t*na = tna, t^cc ^ ^n'a'' ^/ 7^ (^'' «')• 

= Y,(l)nafl{tna){(t)na\^), 
na 

= Y,(t)naf2{tna){(t>na\^), 
na 

Suppose that F is an operator with non-degenerate spectrum; then 

if the observable T with non- degenerate spectrum has a definite value in the state p, then 
it is possible to represent the density matrix of this state as a function of the operator F. 

The operator F can be written in the form 

F — fnPni 
n 

Pn Fni PmPn ^mnPmi tri^P^ 1, ^ ] Pn F. 

n 

The numbers {/„} satisfy the conditions 

fn = fn, fny^fn', ^f n ^ U . 

We obviously have 

F = E/"^' 

n 

Prom 

(F) = J2Pnfn -- 

n 

n 
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we get 

Pn — ^nN- 

In this case density matrix is a projective operator satisfying the condition 

= p. 

It acts as 

where |^) is a vector in Hilbcrt space. 

The average value of an arbitrary variable in this state is equal to 

{A) = {^n\A^n). 

It is so-called PURE state. If the state is not pure it is known as mixed. 

Suppose that every vector in 7i is a square integrable function ^'(x), where a; is a set of 
continuous and discrete variables. Scalar product is defined by the formula 

(^1$) = J dx-^*{x)^x). 

For simplicity we assume that every operator F in Ti. acts as follows . 

(F*)(a;) = J F{x,x')dx'^{x'). 

That is for any operator F there is an integral kernel F{x,x') associated with this operator 

F <^ F{x,x). 

Certainly, we may use (5-function if necessary. 

Now the average value of the variable J-' in the state p is given by equation 

{J-)p — J F{x,x)dx p{x ,x)dx. 

Here the kernel p{x,x') satisfies the conditions 

p*{x,x') = p{x',x), 

J p{x,x)dx — 1, 
V^* e H j '^{x)dxp{x,x)dx'if{x) > 0. 
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3 Composite System and Reduced Density Matrix 

Suppose the variables x are divided into two parts: x — {y, z}. Suppose also that the space 
7i is a direct product of two spaces Tii, H^- 

n = Hit^m. 
So, there is a basis in the space that can be written in the form 

0an(Z/, Z) = fa{y)Vn{z) ■ 

The kernel of operator F in this basis looks like 

F <S=^ F{y,z;y',z) . 
In quantum mechanics it means that the system 5" is a unification of two subsystems Si and 

S = SiLiS2. 

The Hilbert space H corresponds to the system 5* and the spaces Hi and H2 correspond to 
the subsystems Si and 5*2 . 

Now suppose that a physical variable J-'i depends on variables y only. The operator that 
corresponds to J-'i has a kernel 

Fi{y,z;y',z') = Fi{y,y')5{z - z) . 

The average value of Fi in the state p is equal to 

{Fi)p = j F{y,y')dy'pi{y',y)dy, 
where the kernel pi is defined by the formula 

Pi{y,y') = / p{y,z;y',z)dz. 
The operator pi satisfies all the properties of Density Matrix in ^i. Indeed, we have 

pi*iy,y') = piiy',y), 
j Piiy,y)dy = i, 

V^'i e Tioo / ^i{y)dypi{y, y')dy'^i{y') > . 

The operator 

Pi = Tr2Pi+2, 

is called Reduced Density Matrix . Thus, the state of the subsystem is defined by 
reduced density matrix. 

The reduced density matrix for the subsystem ^2 is defined analogously. 

P2 = Tripi+2- 
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Quantum states pi and p2 of subsystems are defined uniquely by the state ^1+2 of the 
composite system. 

Suppose the system S" is in a pure state then a quantum state of the subsystem 5*1 is 
defined by the kernel 

Pi{y,y') = J^{y,z)dz^*{y',z). 
If the function "^{y, z) is the product 

^{y,z) = f{y)w{z), J \w{z)\^dz = l, 
then subsystem Si is a pure state , too 

Pi{y,y') = f{y)ny), j\f{y)?dy = i. 

As it was proved by von Neumann, it is the only case when purity of composite system is 
inherited by its subsystems. 

Let us consider an example of a system in a pure state having subsystems in mixed states. 
Let the wave function of composite system be 

^{y,z) = -^{f{y)w{z)±f{z)w{y)), 

where < f\w >= and < f\f >=< w\w >= 1. The density matrix of the subsystem 5*1 has 
the kernel ^ 

pi{y,y') = ^{f{y)f*{y') +w{y)w*{y')). 

The kernel of the operator pf has the form 

Pi\y,y') = \{f{y)ny')+w{y)w*iy')). 

Therefore, the subsystem Si is in the mixed state. Moreover, its density matrix is propor- 
tional to unity operator. The previous property resolves the perplexities connected with 
Einstein - Podolsky - Rosen paradox. 

4 EPR - paradox 

Anyway, it was Shrodinger who introduced a term " EPR-paradox" . The authors of EPR 
themselves always considered their article as a demonstration of inconsistency of present to 
them quantum mechanics rather than a particular curiosity. 

The main conclusion of the paper [Q "Can Quantum- Mechanical Description of Physical 
Reality Be Considered Complete?" published in 1935 (8 years later then the von Neumann 
book) is the statement: 

..we proved that (1) the quantum mechanical description of reality given by wave func- 
tions not complete or (2) when the operators corresponding to two physical quantities do 
not commute the two quantities cannot have simultaneous reality. Starting then with the 
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assumption that the wave function does give a complete description of the physical reality, 
we arrived at the conclusion that two physical quantities, with noncommuting operators, can 
have simultaneous reality. Thus the negation of (1) leads to negation of only other alternative 
(2). We can thus focused to conclude that the quantum-mechanical description of physical 
reality given by wave function is not complete. 

After von Neumann's works this statement appears obvious. However, in order to clarify 
this point of view completely we must understand what is "the physical reality" in EPR. In 
EPR-paper the physical reality is defined as: 

//, without in any way disturbing a system, we can predict with certainty (i.e., with 
probability equal to unity) the value of physical quantity, then there exists an element of 
physical reality corresponding to this physical quantity. 

Such definition of physical reality is a step back as compared to von Neumann's definition. 
By EPR definition, the state is actual only when at least one observable has an exact value. 
This point of view is incomplete and leads to inconsistency. 

When a subsystem is separated "the loss of observables" results directly from the defi- 
nition of density matrix for the subsystem. "The occurrence" of observables in the chosen 
subsystem when the quantities are measured in another "subsidiary" subsystem can be nat- 
urally explained in the terms of conditional density matrix. 



5 Conditional Density Matrix 

The average value of a variable with the kernel 

F''{x,x') = F,{y,y')u{z)u*{z'), J \u{z)\^dz = 1, 

is equal to 

= P j Fi{y,y)dy' p\y\y)dy, 

where 

P^iy.y) = - / u*{z)dzp{y,z]y',z')u{z')dz' , 
p J 

p ^ j u*{z)dz p{y,z;y,z')u{z')dz'dy. 
Since we can represent p in the form 

p — j P{z, z')dz' p2{z'; z)dz, 

P{z,z') = u{z)u*{z'), 

we see that p is an average value of a variable P of the subsystem 5*2. Operator P is a 
projector (P^ = P). Therefore it is possible to consider the value p as a, probability. 

It is easy to demonstrate that the operator p'^ satisfies all the properties of density matrix. 
So the kernel p'^{y, y') defines some state of the subsystem ^i. What is this state? 
According to the decomposition of (5-function 

5{Z - Z) ^^(t)n{z)(l>n{z), 
n 
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{(I)n{z)} being a basis in the space H2, the reduced density matrix is represented in the form 
of the sum 

pi{y,y') = Y.PnPniy^y')- 



Here 



and 



Pniy^y') = — I (pn{z)dzp{y,z]y',z')(f)n{z')dz' 
Pn = [ (j)„*iz)dz p{y,z;y,z')(pniz')dz'dy 



The numbers p„ satisfy the conditions 

Pn* = Pn, Pn > 0, J^P^ = 

n 

and are connected with a probabihty distribution. 

The basis {4>n} in the space H2 corresponds to some observable G2 of the subsystem 5*2 
with discrete non-degenerate spectrum. It is determined by the kernel 

G2{,Z, Z) = J2 9n<Pn<P*n, 9n = g*n, 9n 7^ 9nl if Til. 
n 

The average value of G2 in the state p2 is equal to 

dyp2{z, z')dz'G{z' , z) = 

= Y.9n / dyp2{z, Z')dz'(l)n{z')(t)n{z') = Y.Pn9n- 
n n 

Thus number pn defines the chance that the observable G2 has the value Qn in the state p2- 
Obviously, the kernel Pn{y, y') in this case defines the state of system 5*1 under condition that 
the value of variable G2 is equal to Qn- Hence it is natural to call operator as Conditional 
Density Matrix (CDM) [|, § 

, _ rr2(Ap) 

It is {conditional) density matrix for the subsystem 5*1 under the condition that the sub- 
system 5*2 is selected in a pure state p2 = A- It is the most important case for quantum 
communication. Conditional density matrix satisfies all the properties of density matrix. 

Conditional density matrix helps to clarify a sense of operations in some finest experi- 
ments. 
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6 Examples: System and Subsystems 



6.1 Parapositronium 

As an example we consider parapositronium, i.e. the system consisting of an electron and 
a positron. The total spin of the system is equal to zero. In this case the nonrelativistic 
approximation is valid and the state vector of the system is represented in the form of the 
product 

*(f;,C7e;fp,(7j,) = $(r;,fp)x((7e,C7p). 

The spin wave function is equal to 

X((7e,(7p) = -^(Xn((^e)X-n((^p) - Xn(f^p)X-n((^e)) • 

Here Xn((j) and X(-n){'^) are the eigenvectors of the operator that projects spin onto the 
vector ft: 

The spin density matrix of the system is determined by the operator with the kernel 
The spin density matrix of the electron is 

In this state the electron is completely unpolarized. 

If an electron passes through polarization filter then the pass probability is independent 
of the filter orientation. The same fact is valid for the positron if its spin state is measured 
independently of the electron. 

Now let us consider quite a different experiment. Namely, the positron passes through the 
polarization filter and the electron polarization is simultaneously measured. The operator 
that projects the positron spin onto the vector m (determined by the filter) is given by the 
kernel 

Now the conditional density matrix of the electron is equal to 
The result of the summation is 

Pe/p (CX, a' ) = X(-m) (CX) X(-m) i^^' ) ■ 
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Thus, if the polarization of the positron is selected with the help of polarizer in the state 
with well defined spin, then the electron appears to be polarized in the opposite direction. 
Of course, this result is in an agreement with the fact that total spin of composite system is 
equal to zero. Nevertheless this natural result can be obtained if positron and electron spins 
are measured simultaneously. In the opposite case, the more simple experiment shows that 
the direction of electron and positron spins are absolutely indefinite. 

A.Eistein said "rajjinert ist der Herr Gott, aber boschaft ist Er nichf . 

6.2 Quantum Photon Teleportation 

In the Innsbruck experiment on a photon state teleportation, the initial state of the system 
is the result of the unification of the pair of photons 1 and 2 being in the antisymmetric 
state x{'^i:'^2) with summary angular momentum equal to zero and the photon 3 being in 
the state x»fi(c"3) (that is, being polarized along the vector m). The joint system state is 
given by the density matrix 

p{a,a) = ^{a)^*{a'), 
where the wave function of the joint system is the product 

= X(0-i,Cr2) Xm(0-3). 

Considering then the photon 2 only (without fixing the states of the photons 1 and 3) we 
find the photon 2 to be completely unpolarized with the density matrix 

p(o-2,cr2) = Tr(i,3) p(o-i, 0-2, 0-3; ai, 4,(13) = ^ ^{(^2 ~ ^'2)- 

However, if the photon 2 is registered when the state of the photons 1 and 3 has been 
determined to be x(^i) '^3) then the state of the photon 2 is given by the conditional density 
matrix 

_ ^^^(1,3) (^1,3 Pl,2,3) 

~ Tr(P,3a2,3) • 

Here Pi^s is the projection operator 

Pi,3 = x(c^i,o-3) X*(o-i,cr3). 
To evaluate the conditional density matrix it is convenient to preliminary find the vectors 

<Pi^l) = X^(^3) X(f^l,^3) 
3 

and 

^(^2) = X{(^1,(^2)- 

1 

The vector 6 equals to 

^(^^2) = -- Xni{(^2) 
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and the conditional density matrix of the photon 2 appears to be equal to 

P2/{1,3} = Xrn{(^2) 

Thus, if the subsystem consisting of the photons 1 and 3 is forced to be in the antisymmetric 
state x(cri,o"3) (with total angular momentum equal to zero) then the photon 2 appears to 
be polarized along the vector m. 



6.3 Entanglement Swapping 

In the recent experiment ^ in installation two pairs of correlated photons are emerged 
simultaneously. The state of the system is described by the wave function 

■q/(a) = ^(0-1,0-2,0-3,0-4) = x(cri,o-2)x(o-3,a4). 

The photons 2 and 3 are selected into antisymmetric state x(o-2,o-3). 
What is the state of pair of photons 1 and 4? 
Conditional density matrix of the pair (1-4) is 

P14/23 — 



Tr{P23Pi234) 

where operator that selects pair (2-3) is defined by 

^23((7,a') = x(0^2,f^3)X*(ff2.O3) 

and density matrix of four photons system is determined by kernel 

Pi234(o-,cr') = ^ (cTi, 0-2, as, (74) **(cri, 0-2, 0-3,^4). 

Direct calculation shows that the pair of the photons (1 and 4) has to be in pure state with 
the wave function 

$(0-1,0-4) = x(cri,o-4)- 
The experiment confirms this prediction. 



6.4 Pairs of Polarized Photons 

Now consider a modification of the Innsbruck experiment. Let there be two pairs of 
photons (1, 2) and (3, 4). Suppose that each pair is in the pure antisymmetric state x- The 
spin part of the density matrix of the total system is given by the equation 

where 

^(o-) = x(o-i,o"2) x(o"3,o^4)- 

If the photons 2 and 4 pass though polarizes, they are polarized along Xrn{cr2) and Xs(o-4) 
then the wave function of the system is transformed into 

^{(^) = Xn{(^l) Xm{(^2) XMs) Xsicr^)- 
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Here n, m and f , s are pairs of mutually orthogonal vectors. 

Now the conditional density matrix of the pair of photons 1 and 3 is 

P(l,3)/(2,4)(^,f^') = ©((^1,^3) ©*(^^1,0^3)- 

The wave function of the pair is the product of wave functions of each photon with definite 
polarization 

Qi<Jl,a3) = Xn(f^l) Xr(cr3). 

We note that initial correlation properties of the system appear only when the photons 
pass though polarizers. Although the wave function of the system seems to be a wave function 
of independent particles the initial correlation exhibits in correlations of polarizations for each 
pair. Pairs of polarized photons appear to be very useful in quantum communication. 

6.5 Quantum Realization of Verman Communication Scheme 

Let us recall the main idea of Vernam communication scheme [0| . In this scheme, Alice 
encrypts her message (a string of bits denoted by the binary number nii) using a randomly 
generated key k. She simply adds each bit of the message with the corresponding bit of the 
key to obtain the scrambled text (s = mi © k, where © denotes the binary addition modulo 
2 without carry). It is then sent to Bob, who decrypts the message by subtracting the key 
{sQk = rui ®kQk = mi). Because the bits of the scrambled text are as random as those of 
the key, they do not contain any information. This cryptosystem is thus provable secure in 
sense of information theory. Actually, today this is the only probably secure cryptosystem! 

Although perfectly secure, the problem with this security is that it is essential that 
Alice and Bob possess a common secret key, which must be at least as long as the message 
itself. They can only use the key for a single encryption. If they used the key more than 
once. Eve could record all of the scrambled messages and start to build up a picture of 
the plain texts and thus also of the key. (If Eve recorded two different messages encrypted 
with the same key, she could add the scrambled text to obtain the sum of the plain texts: 
Si © S2 = mi © A; © m2 © = mi © m2 © A; © = mi © m2, where we used the fact that © 
is commutative.) Furthermore, the key has to be transmitted by some trusted means, such 
as a courier, or through a personal meeting between Alice and Bob. This procedure may be 
complex and expensive, and even may lead to a loophole in the system. 

With the help of pairs of polarized photons we can overcome the shortcomings of the clas- 
sical realization of Vernam scheme. Suppose Alice sends to Bob pairs of polarized photons 
obtained according to the rules described in the previous section. Note that the concrete 
photons' polarizations are set up in Alice's laboratory and Eve does not know them. If 
the polarization of the photon 1 is set up by a random binary number pi and the polar- 
ization of the photon 3 is set up by a number rrii © pi then each photon (when considered 
separately) does not carry any information. However, Bob after obtaining these photons 
can add corresponding binary numbers and get the number rrii containing the information 
{nii ®Pi®Pi = rrii). 

In this scheme, a secret code is created during the process of sending and is transferred 
to Bob together with the information. It makes the usage of the scheme completely secure. 
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7 Conclusion 



Provided that the subsystem 5*2 of composite quantum system 5" = + 5*2 is selected (or 
will be selected) in a pure state P„ the quantum state of subsystem Si is conditional density 
matrix pic/2n- Reduced density matrix pi is connected with conditional density matrices by 
an expansion: 

Pi = ^PnPln/2n] 

here 

The coefficients p„ are probabilities to find subsystem 5*2 in pure states P„. 
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